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Abstract 

We first note it is possible to construct two linear operators defined on two different domains, yet sharing the same 
spectrum using a geometric transform. However, one of these two operators will necessarily have spatially varying, 
matrix valued, coefficients. This mathematical property can be used in the design of metamaterials whereby two 
different domains behave in the same electromagnetic, acoustic, or hydrodynamic way (mimetism). To illustrate 
this property, we describe a feasible invisibility carpet for linear surface liquid waves in a channel. This structured 
metamaterial bends surface waves over a finite interval of Hertz frequencies. 

To cite this article: G. dupont, S. Guenneau, S. Enoch, C. R. Mecanique xxx (2011). 

Resume 

Tapis d'invisibilite dans un canal avec un fluide structure Nous observons en premier lieu qu'il est possible 
de construire deux operateurs lineaires definis sur deux domaines distincts mais qui possedent le meme spectre 
de valeurs propres, par le truchement d'une transformation geometrique. Neanmoins, un des deux operateurs 
aura necessairement des coefficients heterogenes et non scalaires. Cette propriete mathematique pent neanmoins 
etre utilisee dans le design de metamateriaux grace auxquels deux objects distincts presentent les memes ca- 
racteristiques acoustique, optique ou hydro dynamique (mimetisme). Pour illustrer notre propos, nous proposons 
un modele realiste de cylindres rigides judicieusement disposes (un metamateriau appele tapis d'invisibilite) qui 
fonctionne sur une plage de frequences hertziennes pour des vagues de faible amplitude dans un canal. 
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Version frangaise abregee 

Nous considerons un probleme spectral modele qui consiste a trouver les couples de valeurs propres A 
et vecteurs propres associes (j) tels que : 

Ai((/)) = -A(/) = A(/), (1) 

dans un domaine Vti borne dans . Nous nous interessons plus particulierement au cas de conditions de 
Neumann au bord du domaine, en vue d'une application a I'acoustique (voir equations 3-5). II est bien 
connu que la resolvante de cet operateur Ai est compacte dans I'espace de Hilbert H^iVti) (par injection 
compacte de H^{Vti) dans L^(r^i)), et done que le spectre (j{Ai) de I'operateur Ai est un ensemble 
discret de valeurs propres reelles positives tendant vers +00 qui peuvent etre rangees par ordre croissant 
(Gram- Schmidt) . 

II est bon de noter que ce probleme mathematique modelise par exemple la recherche de modes suscep- 
tibles de se propager dans un guide (acoustique, electromagnet ique, hydraulique...) en regime harmonique, 
auquel cas la racine carree de la valeur propre A a la dimension physique d'une frequence par une vitesse. 




a(Ai) = o(A2) o(An) 

Figure 1. If r2i and Q.2 are two bounded domains that can be mapped onto one another, via a change of coordinates 
described by the Jacobian matrix J, two self-adjoint bounded operators Ai = —A and A2 = —T^^^V • T~"'^V respectively 
defined in £(L^(r2i)) and C{L'^{Q2)) have identical spectra cr(Ai) = (7(^2), where T^"*^ is the upper left block of the inverse 
of the symmetric matrix T = J^J/det(J). Our proposal of cloaking is to asymptotically approach the spectrum (7(^2) with 
a sequence of spectra cr(Ary) associated with media of typical heterogeneity size 77, when 77 goes to zero. 



La question que Ton se pose est de savoir si Ton pent construire un autre operateur A2 agissant sur 
un domaine borne Q2 distinct de Qi dont le spectre (j{A2) est identique au precedent. La reponse est 
affirmative dans la mesure ou Ton procede a un changement de variables qui applique le domaine Qi sur 
le domaine 1^2, cf. Figure 1. En effet, les couples de valeurs propres /3 et vecteurs propres associes tjj tels 
que : 

A2W = -T3-3IV • T-^VV^ = (3^ , (2) 

ou les valeurs propres /3 sont reelles positives (T est symetrique), peuvent-etre mises en correspondance 
(une a une) avec les valeurs propres A. 

Ce tour de passe-passe est bien connu des specialistes de I'inversion [1,2] dans le cadre de I'etude 
du probleme de Calderon (qui revient a connaitre les proprietes de I'application Dirichlet-Neumann 
[3]). Neanmoins, nous n'avons trouve nulle part dans la litterature physique un expose mathematique 
elementaire qui prend la mesure de ces implications pour les problemes aux valeurs propres dans les 
resonances de cavites : deux operateurs lineaires definis sur des domaines bornes distincts peuvent 
presenter des spectres identiques, pourvu que Tun au moins ait des coefficients heterogenes anisotropes, ce 
qui renvoie a la question de la reconnaissance de forme d'un tambour a travers sa signature acoustique [4]. 

Les aspects mathematiques sous-jacents depassent le cadre de cette etude, mais nous presentons dans 
la Figure 2 et le Tableau 1 des resultats numeriques qui appuient notre propos : les graphes de gauche 
de la figure 2 et les trois premieres colonnes du tableau 1 de gauche demontrent qu'un canal a vagues 
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Table 1 



Eigenvalues (a) and associated norms of eigenfields (b) : Channel filled with straight boundary filled with homogeneous 
isotropic fluid (flrst column), curved boundary (second coumun) ; Channel with curved boundary fllled with heterogeneous 
anisotropic fluid (third column), Channel with curved boundary filled with structured homogeneous isotropic fiuid (third 
column). We note that the spectrum of the curved channel with transformed fiuid coincides with that of the straight channel 
with homogeneous isotropic fiuid, while the spectrum associated with the structured fiuid is different ; However, the norm 
of the eigenfield is very similar in all three cases. The case of a curved channel filled with homogeneous isotropic fiuid leads 
to different eigenvalues and eigenfields. The fact that we recover the L'^ norm of the eigenfield with the structured fiuid 
is more important for practical applications (the inherent frequency shift can be disregarded if the priority is to bend the 
wavefront of, say, linear surface water waves). 

a bord droit a le meme spectre (aux erreurs numeriques pres) qu'un canal a bord courbe avec un fluide 
transforme (i.e. heterogene anisotrope deduit de la formule (7), voir figure 3). Les graphes de droite de 
la figure 2 et les trois premieres colonnes du tableau 1 de droite demontrent que les normes dans des 
fonctions propres coincident (aux erreurs numeriques pres) pour les canal a vagues a bord droit et a bord 
courbe avec un fluide transforme. 

La version anglaise de Particle est quant a elle dediee a I'etude d'un probleme de diffraction en hy- 
drodynamique qui est le pendant des problemes aux valeurs propres analyses ci-dessus. L'accent est mis 
sur une simulation numerique avec un design de tapis structure qui doit faire I'objet d'une validation 
experiment ale ulterieure dans un canal a houle. II est interessant de noter que les valeurs propres du 



3 



♦ Channel with curved boundary 

* Channel with straight boundary 

• Transformed fluid 
■ Structured fluid 
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♦ Channel with curved boundary 

* Channel with straight boundary 

• Transformed fluid 
■ Structured fluid 



Figure 2. Graph of eigenvalues (a) and associated norms of eigenfields (b) : Panel (a) The red stars for the straight 
channel are superimposed with the blue circles for the curved channel with transformed fluid. The green diamonds for the 
curved channel with isotropic fluid and the black squares for the stuctured fluid are both far away from the red stars. 
Panel (b) : The red stars for the straight channel are superimposed with the blue circles for the curved channel with 
transformed fluid. For most points, the green diamonds (curved channel with isotropic fluid) are further away from red stars 
and blue circles than the black squares (stuctured fluid) : The eigenfleld shares similar features within the structured and 
the transformed fluids. 

probleme spectral associees au tapis structure different de celles du meta-fluide heterogene anisotrope ob- 
tenu par transformation geometrique, cf. la troisieme colonne du tableau 1 (gauche), alors que les fonctions 
propres correspondantes coincident presque, cf. la troisieme colonne du tableau l(droite). Ces resultats 
sont consistants avec les cartes de champ associees a ces deux cas, cf. figures 3, 4 et 5. L'ensemble de 
ces resultats numeriques suggere que le design de cloques structurees dans des canaux a vagues necessite 
un travail ulterieur d'optimisation (de type probleme inverse) sur des problemes spectraux afin de mieux 
cerner leurs proprietes intrinseques : les interations fluides-structures conduisent au mimetisme a certaines 
frequences (au sens ou les vagues dans un canal droit sont similaires a celles d'un canal courbe), mais aussi 
a des resonances propres du canal a d'autres frequences (auquel cas le cloaking est caduc). Ce dernier 
point s'apparente aux modes de resonance des capes d'invisibilite decouverts par I'equipe de Greenleaf 
dans le contexte de la mecanique quantique [2], qui ont aussi leur pendant en acoustique [5]. 



1. Setup of the hydrodynamic problem 

The transformation based solutions to the Maxwell equations in curvilinear coordinate systems reported 
by Pendry et al. in [6] bend electromagnetic waves around arbitrarily sized and shaped surfaces (see also [7] 
for a conformal optics approach). The electromagnetic carpet is a metamaterial which maps a concealment 
region onto a surrounding surface: as a result of the coordinate transformation the permittivity and 
permeability are strongly heterogeneous and anisotropic within what physicists call an invisibility carpet 
[8] , yet fullling impedance matching with the surrounding vacuum. The carpet thus neither scatters waves 
nor induces a shadow in the reflecting field. 

In the present paper, we build upon the recent proposal by Li and Pendry [8] to map a curved surface 
onto a flat surface in order to control the wave front of an electromagnetic wave scattered by a bump 
located on a flat mirror so that if we now dress this bump with a heterogeneous anisotropic material, 
the wave seems to be reflected by a flat mirror, thereby making the bump invisible [8,9]. We actually 
design a structured material in order to mimick the prerequisite material properties deduced from a 
geometric transformation. We focus here on linear surface water waves propagating within a channel, but 
we emphasize that our approach is generic and works for any wave governed by a Helmholtz equation 
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subject to Neumann boundary conditions e.g. anti-plane shear waves in an elastic material with cracks, 
pressure water waves in a fluid with rigid inclusions, transverse electric waves in a dielectric medium with 
infinite conducting inclusions. 

Let denote the region of a channel occupied by a fluid. The conservation of momentum leads to the 
Navier- Stokes equations: 



where u denotes the velocity field, t the time variable, p the density of the fiuid /i its viscosity, and g the 
gravity. 

If we assume that the fiuid is incompressible and irrotational, we know that u derives from a potential 
which under the hypothesis of small perturbations of the free interface separating the fiuid with ambient 
atmosphere, leads to the Helmholtz equation: 



with the spectral parameter related to the frequency of the wave uj via the dispersion relation: 



Here, h is the depth of water in the channel and a the surface tension at the free surface. 
The linearized problem (4-5) allows for straightforward analogies between transverse electromagnetic and 
acoustic waves propagating in structured cylindrical domains, see [10] for the design of an invisibility 
cloak for surface liquid waves, experimentally shown to work between 10 and 15 Hertz (broadband). 



2. Design of a heterogeneous anisotropic fluid 

As we already announced in the French abridged version, our aim here is to approximate the spectrum 
of the Laplace operator defined on a bounded region of a certain shape with a perturbed Laplace operator 
defined on another bounded region of a different shape. In both cases, we assume Neumann boundary 
conditions, so that the resolvents of both operators are compact and their spectra consist only of a 
countable set of discrete eigenvalues with a single accumulation point (0 or infinity depending upon 
whether we look at the operator or its inverse) [11]. This allows for a one-to-one correspondence between 
the spectra of the Laplace operator Ai and perturbed Laplace operators Arj associated with the structured 
fiuid. The underlying asymptotic mechanism is that one wants to approximate each eigenvalue of by an 
eigenvalue of A^, in the limit when rj goes to zero: The smaller 77, the larger the number of rigid cylinders 
(of order r]~^) of decreasing diameter (^ 7^), the finer the approximation (in the homogenization limit) 
of the transformed fiuid obtained by mapping the first region on the seconde one. This means that the 
sequence of spectra cr{Arj) of operators Arj should tend (pointwise) to the spectrum (7{Ai) of the operator 
Ai when 77 tends to zero. As a result, if the first region Oi is filled with an isotropic homogeneous medium 
(say a fiuid), the second one Cl2 (associated with the domain of an operator A2) is now filled with the same 
fiuid, however with a collection of small rigid cylinders approximating an anisotropic heterogeneous fiuid: 
(j{Ai) = limrj^Q a {Arj) = (j{A2). We now want to numerically validate this conjecture: one can design 
a meta- fluid by structuring the second region with rigid cylinders, in which case it can be simply filled 
with an ordinary fiuid. However, such an asymptotic approach can only work to certain extent (within 
the framework of effective medium theory, hence for small enough frequencies). 




(3) 



(4) 




(5) 
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2.1. Geometric transform 



Let us first introduce a simple geometric transform mapping the first region to the second one. The 
bottom line is the bold proposal by Li and Pendry to conceal an object that is placed under a curved 
reflecting surface by imitating the reflection of a flat surface [8] in the context of electromagnetic waves 
in open space. In the present case, the domain is bounded and the geometric transform reads as follows: 
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(6) 



and a = {y2- yi)/yi- 

The metric tensor associated with the transformed coordinates takes the following form (and its effect on 
the Cartesian metric is shown in figure 3): 




(7) 
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Figure 3. Left: Metrics associated with the Cartesian coordinate system (original domain, leftmost panel) and the transformed 
coordinate system (invisibility carpet, right panel) mapped onto one another via the transformation matrix T (note that 
the right angles are not preserved i.e. the transformation is not conformal); Right: Numerical simulations at frequency 
u = 1.99Hz; (a) Field inside a straight channel filled with a homogeneous isotropic fluid; (b) Field inside a curved channel 
fllled with a homogeneous isotropic fluid; (c) Field inside a curved channel fllled with a heterogeneous anisotropic fluid 
described by formula (7). The color scale is in arbitrary units. The strong similarity between flelds in (a) and (c) is noted. 

It is interesting to look at the expression of the eigenvalues of as these are the relevant quantities to 
design a structured channel: 




-4a2 + 1 + a2 



dy 
dx' 



(8) 



> 



We note that Ai and A^, i = 2,3, are strictly positive functions as obviously 1 + + 
/— + ^1 + + ^^^^ a > 0. This establishes that is not a singular matrix for 
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a two-dimensional carpet, which is a big advantage over two-dimensional cloaks obtained by blowing up 
a point onto a disc [8]: the transformation matrix is then singular at the cloak's inner boundary (one 
eigenvalue goes to infinity, while the other two go to zero [3]). 

2.2. Structured fluid 

Let us now mimic the heterogeneous anisotropic fiuid using an effective medium approach whereby 
an assembly of rigid cylinders judiciously located is now fixed to the bottom of the channel. It is clear 
that such a design will only work to certain extent and moreover will be constrained by the working 
eigenfrequency (the larger the eigenvalue of the operator, the larger the discrepancy between the ideal 
and approximated cases). In Figures 4 and 5, we show some representative fields corresponding to given 
eigenfrequencies in the range 1.72Hz < v < 2.33Hz for a curved channel with a carpet (Figure 4) and 
without a carpet (Figure 5). We emphasize that the wavefront of the fields is nearly fiat in Figure 4. We 
report in table 1 the norm of these eigenfields and compare them to the benchmark of a straight channel 
and a curved channel filled with a transformed fiuid. These numerical results clearly show the positive 
effect of the structured carpet. However, some care needs be taken when commenting these results, as 
shown by the discrepancy between the eigenvalues for spectral problems set in the straight channel and 
the curved channel filled with the structured fluid: cloaking is only achieved for the control of the fields's 
wavefront, not for the eigenfrequencies. 




Figure 4. Eigenfields for 1.72 Hz < v < 2.23 Hz in a curved channel with the structured carpet. The flat wavefronts of all 
eigenflelds is noted. 
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Figure 5. Eigenfields for 1.72 Hz < v < 2.23 Hz in the same curved channel as in figure 4 but without the structured carpet. 
The disturbed wavefronts for the eigenfields is noted (except for the first two leftmost eigenfields). 

3. Conclusions 

In this note, we have reported some prehminary results on a structured invisibihty carpet for the 
control of linear surface water waves in a channel. Unlike for the structured invisibility cloak some of us 
designed earlier for linear surface water waves [10] (which avoids any backscattering of an incident wave), 
the carpet mimics the backscattering of a flat boundary (i.e. it flattens the wavefront of backscattered 
waves). The numerical illustrations demonstrate the high potential for a practical realization of a meta- 
fluid working over a large bandwidth. We hope this analysis will foster experimental efforts towards a 
new generation of dykes without overtopping phenomena. Similar ideas could be implemented in the 
design of structured fluids for an enhanced control of pressure waves [12]. It should be finally pointed out 
that recent theoretical and experimental work drawing analogies between water waves and cosmological 
physics [13,14] suggests new avenues for structured meta- fluids in the non-linear regime. 
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